求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

Permutation trinomials over Fq3

We consider four classes of polynomials over the fields ${\mathbb{F}}_{q^3}$, $q=p^h$, $p>3$, $f_1(x)=x^{q^2+q-1}+A{x}^{q^2-q+1}+Bx$, $f_2(x)=x^{q^2+q-1}+A{x}^{q^3-q^2+q}+Bx$, $f_3(x)=x^{q^2+q-1}+A{x}^{q^2}-Bx$, $f_4(x)=x^{q^2+q-1}+A{x}^q-Bx$, where $A,B\in {\mathbb{F}}_q$. We find sufficient conditions on the pairs $(A,B)$ for which these polynomials permute ${\mathbb{F}}_{q^3}$ and we give lower bounds on the number of such pairs.     

基于LBM方法的风力发电机尾流结构仿真

针对传统CFD数值计算方法难以实现风力机动态旋转及其旋转状态下的流固耦合计算,本文结合格子玻尔兹曼(LBM)方法易于处理动态复杂边界的特点及大涡模拟(LES)方法在非稳态涡流结构捕捉上的优势,采用LBM-LES联合方法进行三维风力发电机整机气动性能及尾流结构仿真研究,同时采用尺度自适应方法对尾涡结构进行跟踪和精细化计算。针对NREL PhaseⅥ型试验机进行模拟,得到了与实验结果吻合的流动形态及尾流结构演变规律,分析了尾流区速度演变规律并对比了不同亚格子湍流模型对计算结果的影响.     

一种超导重力仪磁悬浮力的等效计算

为实现超导重力仪磁悬浮力的精确计算,以GWR型超导重力仪为模型基础,采用有限元的思想,将超导球表面电流理想化为多个等高共轴电流环,计算出各个电流环与超导线圈的作用力,求和得到线圈与超导球间的磁悬浮力。利用MATLAB完成计算程序实现,通过改变下线圈电流和上、下线圈电流比,获得满足一定条件的磁悬浮力及其梯度。选取合适的模型参数,计算出线圈对质量为m=4.069 g超导球的磁悬浮力大小为:Ftotal=3.988×10^-2N,磁悬浮力梯度为:-9.699×10^-3N/m,此时悬浮力梯度合适,满足系统稳定性和灵敏度的要求。